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Abstract 

The leading electronic interactions governing the local physics of the CE phase of half-doped 
manganites are extracted from correlated ab initio calculations performed on an embedded cluster. 
The electronic structure of the low-energy states is dominated by double exchange configurations 
and 0-2p a to Mn-3d charge transfer configurations. The model spectra of both a purely magnetic 
non-symmetric Heisenberg Hamiltonian involving a magnetic oxygen and two non-symmetric dou- 
ble exchange models are compared to the ab initio one. While a satisfactory agreement between 
the Heisenberg spectrum and the calculated one is obtained, the best description is provided by 
a double exchange model involving excited non-Hund atomic states. This refined model not only 
perfectly reproduces the spectrum of the embedded cluster in the crystal geometry, but also gives 
a full description of the local double- well potential energy curve of the ground state (resulting from 
the interaction of the charge localized electronic configurations) and the local potential energy 
curves of all excited states ruled by the double exchange mechanism. 

PACS numbers: 75.47.Lx, 71.10.-w, 75.30.Mb, 75.30.Et 
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I. INTRODUCTION 



Although numerous theoretical and experimental studies have been devoted to the un- 
derstanding of the various phase transitions occuring in doped manganites,- their electronic 
structure near half doping is still not fully understood.- Lao.sCao.sMnOa is antiferromag- 
netic below the Neel temperature T/v and paramagnetic for T > Tco, where CO stands for 
charge-ordered. For intermediate temperatures T/v < T < Tco t wo different electronic stuc- 
tures have been proposed. In the first one, a charge order prevails and the magnetic phase 
exhibits a CE-type ordering of the magnetic moments.-! 4 - In the second electronic structure,- 
the holes are trapped in pairs of equivalent Mn ions of intermediate valence state Mn 3,5+ 
resulting from the double exchange (DE) interaction between the two limit Mn 3+ 2_ Mn 4+ 
and Mn 4+ 2_ Mn 3+ electronic configurations. The latter electronic structure is referred as 
a Zener polaron order in reference to Zener who first proposed the DE mechanism.- 

Ab initio calculations confirm the existence^ of a charge-ordered phase (partial localiza- 
tion of the holes on one type of Mn ions) in the crystallographic structure of Radaelli et al— 
while a polaronic order-^^ (delocalized holes between pairs of Mn ions) is obtained for the 
crystallographic structure of Daoud Aladine et al— It is worth to note that in both cases cor- 
related ab initio calculations point out that there is a rather strong O to Mn charge transfer, 
resulting in a partial localization of the holes on the bridging oxygens. UHF and DFT peri- 
odic calculations 1 ^^ 3 - on Lao.sCao.sMnOs even lead to a dominantly Mn 3+ 0~Mn 3+ ground 
state charge distribution. 

The impact of this partial localization of the holes on the oxygen sites has recently been 
studied in the peculiar case of a Zener polaron i.e. using the crystal structure of Daoud 
Aladine et a/.- The two Mn ions in the Zener polaron have an almost identical coordination, 
and hence, several symmetric models have been derived and their solutions contrasted with 
ab initio results.-^ It has been shown that i) a truncated-Hubbard in which all the dominant 
configurations are treated variationally can be mapped both on the Heisenberg and the 
usual DE model owing to the mixing of the bridging oxygen and Mn ions orbitals.- ii) 
the energies of the eight lowest states of the Heisenberg and the usual DE models are 
analytically identical and in qualitative agreement with those of the ab initio spectrum.-^ 
iii) a significant improvement is observed when the non-Hund excited atomic states are 
explicitely introduced in the modelization.- The remarkable accuracy of the resulting model 
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spectrum finally shows that the Zener polaron physics is actually ruled by a refined DE 
mechanism where non-Hund atomic states play a non negligible role. 

The present paper focuses on the local electronic structure of the charge-ordered crystallo- 
graphic structure of Lao.sCao.sMnOa. For this purpose we have studied an embedded cluster 
containing two corner-sharing MnOg octahedra using the crystal structure determined by 
Radealli et at.— The two Mn ions in the cluster have different coordination spheres. One site 
is Jahn- Teller distorted, while the other has six (almost) equivalent, rather short Mn-0 dis- 
tances. The two sites are commonly interpreted as occupied by Mn 3+ (3<i 4 ) and Mn 4+ (3c/ 3 ), 
respectively. According to the position of the holes, different electronic orders are expected: 

• If the holes are localized on the Mn sites, the resulting seven unpaired electron sys- 
tem is of mixed valence nature. The DE mechanism induces a resonance between 
Mn 3+ Q 2- Mn 4+ an( j Mn 4+ 2 "Mn 3+ which in that peculiar case will result in a charge 
ordering of the electronic structure i.e. a partial localization of the e s -like extra elec- 
tron (or hole) on one Mn ion Mn 3,5+5 2 ~Mn 3,5_<5 . The t2 S -like electrons are unpaired 
and strongly localized on each ion. 

• If the holes are localized on the bridging oxygen, the corresponding charge distribution 
suggests a dominant purely magnetic local order (Mn 3+ 0~Mn 3+ ) in which the cluster 
would be a ferrimagnetic entity involving a magnetic oxygen and therefore nine un- 
paired electrons. The model Hamiltonian that provides a relevant description of such 
a local electronic order is a Heisenberg Hamiltonian involving non-symmetric exchange 
integrals between each Mn ion and the magnetic oxygen. 

In order to decide which Hamiltonian is most appropriate to describe the local electronic 
structure of this non-symmetric two-center cluster, we compare the different model spectra 
with the ab initio one. The here-used procedure consists in extracting the local effective 
interactions of the model Hamiltonians from calculations performed with the exact electronic 
Hamiltonian. For this purpose, the low energy spectrum of an embedded cluster of the 
material is studied by means of correlated ab initio calculations. The subsequent extraction 
of the model parameters makes use of the effective Hamiltonian theory of Bloch.— ^ The 
comparison of the resulting model spectra to the ab initio one provides rational arguments 
to determine the most appropriate model Hamiltonian. 
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Section II is devoted to the presentation of the embedded cluster computational proce- 
dure. The Heisenberg and DE models as well as their analytical solutions are derived in 
section III. Finally, the comparison of the model spectra with the ab initio one is discussed 
in the fourth section. 



II. DESCRIPTION OF THE CLUSTER AND COMPUTATIONAL PROCEDURE 

The geometry of the M^On cluster is derived from the 20K modulated structure of 
Lao.sCao.sMnOa reported in Table I of Ref. |4j. Two different types of Mn sites can be 
distinguished in the material, the mean radius of their coordination spheres are respectively 
1.971 A and 1.916 A. The distances between a bridging oxygen and the two distinct Mn 
sites are 1.915 A for Mn a -O c and 2.068 A for Mnb-O c . The cluster is represented in Fig. [TJ 
and is constituted of two different Mn ions in a distorted octahedral environment bridged by 
O c , and ten external oxygen atoms. Axis are chosen such that the origin of the coordinates 
is at half distance between the two Mn sites, z is the intermetallic axis and the bridging O 
atom (O c ) lies on the xz-plane. 




FIG. 1: (color online) Mn20n cluster embedded in Mn and (La,Ca) total ion potentials (TIP's). 
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To reproduce the crystal environment, the cluster is embedded in a set of optimized 
point charges which accurately reproduces the Madelung potential of the crystal. Total Ion 
Potentials (TIP) have been used to represent the Mn, La and Ca cations coordinated to the 
oxygen atoms of the cluster (see Fig. [T]). Additional information can be found in Ref. ?], 
where it is shown that this embedding provides reliable results on the ground state of the 
system. 

The electronic structure of the cluster has been studied using the complete active space 
second-order perturbation theory (CASPT2) methocr^^^l implemented in the MOLCAS 
package.— The orbitals are expanded with extended basis sets of ANO type:— & (5s, 4p, 
3d, If) for Mn atoms, (4s, 3p, Id) for the bridging O and (4s, 3p) for external O. The 
CASPT2 method, which has successfully been applied to the study of many spectroscopic 
problems, introduces both non-dynamical and dynamical correlation effects. In a first step, 
a multireference wave function is optimized by the complete active space self-consistent field 
(CASSCF) procedure. This A r -electron wave function includes the main non-dynamical 
electron correlation. The here- considered active space contains all electronic configurations 
obtained by distributing nine electrons in the ten 3c? orbitals of the two Mn ions and the 2p a 
orbital essentially localized on the bridging oxygen. Accurate relative energies of the ground 
state and lowest excited electronic states are computed by introducing the dynamical electron 
correlation in a subsequent second-order perturbation treatment applied to the zeroth-order 
CASSCF wave function. 

Using the ab initio spectrum and the effective Hamiltonian theory of Bloch,— ^ it is 
possible to extract the main electronic interactions of the different models. Since the number 
of electronic interactions is smaller than the number of electronic states described by the 
model space, the value of the effective interactions is fitted to optimally reproduce the ab 
initio spectrum. The mean error per state can be evaluated with the following expression: 

100 x £. \Ef - E?\ 

N x AE AI 1 ' 

where N is the number of states described by the model Hamiltonian, E AI and E M are 
respectively the ab initio and model energies of state i, and AE AI the ab initio spectrum 
width. 

Both the comparison of the so-obtained model spectra to the ab initio one and the size of 
the error per state give rational arguments to discriminate between different modelizations. 
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III. DESCRIPTION OF THE DIFFERENT MODELS 



A. The Heisenberg model 

If the holes are essentially localized on the oxygens O c then the wave functions of the 
lower energy states are dominated by electronic configurations with nine unpaired electrons. 
This is compatible with a Mn^ + -0^-Mn^ + electronic struture and the appropriate model 
Hamiltonian is the Heisenberg Hamiltonian. In a non-symmetric cluster, its specific form as 
a function of the exchange integrals between the three magnetic centers is: 

H Hels = — JiSjVfria ' So c ~~ ^^A/nt " So c — J2^Mn a ' Slfn b , (2) 

where J± is the magnetic interaction between Mn a and O c , J3 the magnetic interaction 
between Mn& and O c , and J2 the magnetic interaction between the two Mn centers. The 
difference between the J\ and J3 exchange integral values reflects not only the geometrical 
asymmetry (the Mn a -O c distance is shorter than the Mnt,-O c distance) but also any possible 
asymmetry in the electronic structure. 

The model space of this Hamiltonian is spanned by the products of the magnetic sites 
ground states. In the considered case, each Mn atom has four unpaired electrons coupled 
to a quintet ground state. The bridging oxygen O c has one unpaired electron in a local 
2po-like orbital leading to a doublet ground state. An analytical solution of the energies 
of the eigenstates of the Heisenberg model space has been derived for this simple case. Its 
expression is: 

E H (S,±) = ±^ [S% ax (S% ax + 1) - S(S + 1)]** + fP (s + fj +J2[S(S+l)-SZ ax (S% ax -l)} 

where x = J± — J3 and (3 = J\ + J3 — 2J 2 . S^ ax is the maximal value that S can take in 
the Heisenberg model. S^ ax = 9/2 in the here-considered cluster. The mean energy of the 
octet states is taken as the zero of energy. 

B. The usual non-symmetric double exchange model 

With the holes localized on the Mn ions, the A-electron wave functions of the lower states 
are dominated by spin configurations presenting seven unpaired electrons, i.e. a Mn^ 5+<5 - 
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O^-Mn^- 5 - 5 (closed shell oxygen) electronic structure. The cluster is of mixed valence 
nature and the DE model should accurately reproduce its physics. 

In a charge-ordered phase, i.e. with 5^0, the electronic structure of the cluster is 
such that the extra electron (or hole) is partially localized on one of the two Mn ions. 
The nuclear relaxation resulting from this partial localization results in an enlargement of 
the coordination sphere (an increase of the average Mn-0 distance) of the Mn ion bearing 
the extra electron and a contraction of the coordination sphere of the Mn ion bearing the 
hole. Piepho, Krausz and Schatz 2 ^ proposed a vibronic coupling model to describe this 
phenomenon in mixed-valence compounds. The relevant dimensionless coordinate in this 
so-called PKS model is q, which describes the nuclear relaxation during the electron transfer 
process. 

In the PKS model q corresponds to the antisymmetric combination of the local stretching 
mode of the coordination spheres. The here-considered coordinate q also includes the Jahn 
Teller distorsion, since the cluster is studied at the experimental geometry which is Jahn- 
Teller distorted. 

Let us call $ a and $6 the electron localized diabatic states corresponding to functions in 
which the extra electron is localized on Mn a and Mn;, respectively. In the usual DE model 
these functions are products of the Mn ions ground state wave functions, i.e. in the m s =7/2 
subspace $ a = Q\-Qb and = Q a -Ql where Q l is a quintet state and Q stands for a quartet 
state. Their energetic dependence to the coordinate q is usually assumed to be quadratic, 
so that the potential energy curves of $ a and obey the equations E a = A(q 2 + 2q) and 
Ef>=A(g 2 — 2q) respectively, where A is the curvature of the diabatic spin states potential 
energy curves as functions of q. The PKS Hamiltonian expressed in the m s =7 /2 subspace 
is : 

A(g 2 + 2q) t \ 
t A(g 2 - 2q) ) 

with t the hopping integral of the extra electron between the two metal centered orbitals 
(e p -like in the present case). Its diagonalization generates the two adiabatic octet states 
0+ = N(® a + $ 6 ) and CL = iV($ a - $ 6 ) of energy 

E PKS (7/2, ±) = Aq 2 ± ^(Ag) 2 + t 2 (3) 

In a non-symmetric and therefore charge-ordered system the potential energy curve of the 
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ground state 0+ as a function of the coordinate q presents a double well. The corresponding 
minima are associated to the two non-symmetric geometries and the Mn 3 ' 5 ~ 5 -0 2_ -Mn 3,5+5 
or Mn 3 5+5 -0 2 ~-Mn 3 ' 5-5 charge-ordered electronic configurations. 

The coupling between the lower m s components of the local ground states $ a and 
generates states of lower spin multiplicity, namely two sextet S- and S + , two quartet Q_ 
and Q + , and two doublet D_ and D + states. In the general case of a non-symmetric 
homonuclear bimetallic complex, the eigenenergies of the model are analytically known 2 ^^ 
and are given by the expression: 



E DE (S, ±) = ±y 4(Ag) 2 + (s + ' + J[S(S + 1) - S™{S™ + 1)] (4) 

where S®^ x is the maximal value that S can take in the DE model, S^ x = 7/2 in the con- 
sidered cluster. B = -^J — r is the leading interaction of the model, t favors a ferromagnetic 
order, in other words, it accounts for the appearance of a high-spin ground state. J is the 
overall exchange integral which describes the electronic circulation in the other open shells, 
here the ^Tike orbitals. This Heisenberg contribution is generally antiferromagnetic and 
has been introduced by Girerd and Papaefthhymiou.— >2I It causes a significant stabilization 
of the mean energies of the low and intermediate spin states with respect to the mean energy 
of the highest spin states. The zero of energy is placed at the mean energy of the highest 
spin states. 

The values of t, J and the product Aq are determined by minimizing e, defined in Eq. [TJ 
The expression of the variable q m in at the minimum of the ground state energy can be 
obtained from the zero of the first derivative of the energy: 

2Ag 

qmm v/4(Ag) 2 + 1 2 1 j 

Knowing the product Aq and the value of q at the minimum, it is possible to determine the 
value of the curvature A and to draw the full potential energy curves as a function of q. 

The hopping parameter can also be extracted directly from the energies and wave func- 
tions of the two octet states without fitting the ab initio spectrum to a model Hamiltonian. 
To validate our numerical fitting procedure, we compare the t-values obtained in both pro- 
cedures. The symmetrically orthonormalized projections of the adiabatic octet state wave 
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functions onto the model space are: 



cos0|$ a ) + sm0|$&) , 



*2> 



sin(j)\& a ) + cos0|$ fe ). 



(6) 



where j—<fr characterizes the charge localization. From these expressions, the des Cloizeaux— 
formalism allows us to extract the hopping integral as follows: 



t = sin(f)cos<p[E{0, -) - E(0, +)] 



(7) 



C. Explicit treatment of the non-Hund excited atomic spin states 

The explicit consideration of the excited atomic non-Hund states in the DE model was 
introduced by Anderson and Hasegawar^ in the description of a symmetric (i.e. fully de- 
localized) cluster. In comparison to the usual DE model, the Anderson-Hasegawa model 
does not only consider the products of local ground states. Its model space is extended to 
products of an atomic ground state on one Mn ion and a single excited non-Hund atomic 
state on the other Mn ion. The so-obtained functions interact with the model space of the 
usual DE model through a term proportional to the hopping integral t. The main contri- 
bution of a variational treatment of these functions is antiferromagnetic. This means that 
the energies of the highest spin states are not affected, while the intermediate and low spin 
states are stabilized by the non-Hund state contribution. As already shown in symmetric 
dimers,- 1 ^^ 1 ^ 1 ^ the best results are obtained using this refined DE model in combination 
with the antiferromagnetic contribution of Girerd-Papaefthymiou through a magnetic ex- 
change integral J. Assuming that the potential energy curves of the electronic states formed 
by a product of a non-Hund state and an atomic ground state are parallel to those of the 
products of the atomic ground states, it is possible to generalize this refined model to the 
study of non-symmetric clusters. Calling Sh the relative energy of the non-Hund state, one 
obtains the following energy expression of the DE states: 



e nh(Si ±) 



2 



1 



5 2 H + lQ(Aq) 2 + At 2 ±A6 H JA(Aq) 2 + B 2 I S+- 




) 



2 



+ J[S(S+1)-S™(S™ + 1)] 



(8) 
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Again, the electronic interactions t, 5h, J and the product Aq are numerically optimized 
such that the model spectrum optimally fits the ab initio one. Since the contribution of 
the non-Hund states is zero in the highest spin states and therefore in the ground state, the 
extraction of q m i n (and then A) is performed using Eq. 

IV. CONTRASTING THE MODELS TO AB INITIO RESULTS 
A. Comparison of the spectra 

Taking the CASSCF wave functions as references, we determined the second-order per- 
turbation theory MS-CASPT2 spectrum for the low-lying states of the cluster. The ordering 
of the states is the following: 

E{OJ) < E{SJ) < E{QJ) < E(D_) < E{D + ) < E{Q + ) < E(S+) < E(0+), (9) 

which is compatible with all the here- considered models. The electronic interactions of the 
different models optimized in order to reproduce at best the ab initio spectrum are given 
in Table [J and the corresponding model spectra are represented in Figure [2] together with 
the MS-CASPT2 spectrum. Using Eq. UJ we validate the numerical fitting procedure. 
Unfortunately, the comparison cannot be done at the MS-CASPT2 level since this method 
does not provide the necessary wave function coefficients.— Instead, we take the spectrum 
provided by the zeroth-order wave function and extract t both by fitting the usual DE 
Hamiltonian and by applying Eq. [7J The respective values of 0.9043 eV and 0.9108 eV 
compare very well. 

From the comparison of the different spectra, it appears that the qualitative features of 
the ab initio spectrum are reproduced in all model spectra, legitimating the consideration 
of the three different models. Notice that this is not the case for all Hamiltonians: The 
mechanism proposed by Zener,- which only considers the hopping integral and does not 
introduce any antiferromagnetic contribution, cannot reproduce the ab initio spectrum.- 

Concerning the extracted electronic interactions and the mean error per state of the 
different models, several comments can be made. 

• The mean error per state of the usual DE and Heisenberg models are similar but not 
strictly identical as in the symmetric cluster.-^ 
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• As already observed for symmetric DE compounds, the explicit consideration of the 
non-Hund atomic states improves quantitatively the reproduction of the ab initio spec- 
trum and therefore the modelization. An error of only 0.14% is obtained for the refined 
DE model. 

• The non-symmetry of the cluster (which is responsible for the charge ordering) is 
reflected in the difference of the J\ and J3 exchange integral values of the Heisenberg 
model. The very large value of J3 shows the covalent character of the interaction 
between Mn& and the bridging O c oxygen. Actually, the value is so large that it even 
questions the validity of a purely magnetic description of the cluster. 

• The values of the electronic interactions t and J extracted from the two DE mod- 
els are qualitatively different. In particular, the value of the exchange integral J 
of the usual DE model is twice as large as the one extracted from the refined DE 
model. A large difference between the two extracted values has already been ob- 
served in other studies.-&22>2i The participation of non-Hund states in the ab initio 
wave function enhances the covalent character of the interaction between the Mn sites. 
Their explicit inclusion in the model space results in an antiferromagnetic contribu- 
tion, i.e. it stabilizes the low and intermediate spin states. In order to reproduce this 
antiferromagnetic contribution, the usual DE model overestimates the value of the 
antiferromagnetic exchange integral J and underestimates the value of t, i.e. of the 
ferromagnetic contribution. 

B. Role of the Non-Hund states on the spin multiplicity of the ground state 

As already shown in the peculiar case of a symmetric cluster,-^ both the exchange integral 
J which describes the electronic circulation in the i2g-like orbitals and the non-Hund states 
stabilize the low and intermediate spin states with respect to the highest spin one. If these 
antiferromagnetic contributions are large enough they may even generate an intermediate 
spin ground state.— ^ We have seen in the previous subsection that the model spectra present 
the same qualitative features as the ab initio spectrum at the experimental geometry. 

As shown in Sec. IHI B\ the curvature A of the potential energy curves of the diabatic 
spin states can be extracted from the calculated spectrum. From the value of A and the 
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TABLE I: Effective parameters of the refined DE (non-Hund), usual DE (Zener, Girerd Pa- 
paefthymiou, ZGP), and Heisenberg model extracted from the MS-CASPT2 spectrum. The mag- 
netic exchange (J) and hopping (t) integrals are in eV. q m in is the value of the effective parameter 
q at the minimum of the ground state potential energy curve, A (in eV) is the curvature of the 
potential energy curves of the diabatic (i.e. before interaction) left and right states, e is the mean 
error of the fit of the model Hamiltonian spectrum to the ab initio spectrum. 



non-Hund ZGP Heisenberg 



Jl 






0.568 


■h 






0.939 


■h 






0.082 


J 


0.0405 


0.0847 




5h 


4.993 






A 


0.7177 


0.6935 




[rain 


±0.38 


±0.32 




t 


1.3292 


1.3081 




C 


0.14 


1.014 


1.010 



optimized effective electronic interactions, it becomes possible to determine the full local 
potential energy curves of the adiabatic states as functions of the coordinate q for the two 
DE models. Figures [3] and @] depict the curves of the eight spin states for the usual and 
the refined DE model, respectively. It is striking to see that while the octet state is the 
ground state at the experimental geometry in the usual DE model, the potential energy 
curves show more stable minima for low and intermediate spin states at larger values of q. 
On the contrary, the octet potential energy curve is the lowest state over the whole interval 
for the refined DE model. 

This qualitative disagreement between the usual DE model and the ab initio calculation 
is due to the overestimation of the exchange integral J. In absence of the hopping integral t 
(i.e. for the diabatic states), the potential wells associated with the left and right localized 
holes are such that the ordering of the eight lowest electronic states is governed by the 
Heisenberg contribution. This antiferromagnetic interaction induces an antiferromagnetic 
order and makes the doublet state the ground state. From Eq. 4 and 7 it is easily seen that 
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■1,5 



Q+ 



FIG. 2: Comparison of the MS-CASPT2 spectrum of the embedded cluster (column 1) with the 
outcomes of the refined DE model (column 2), the usual double exchange model (colum 3) and the 
Heisenberg model (column 4). 

the coupling between the left and right diabatic states through the hopping integral is more 
effective in the highest spin state. Hence, this interaction, which is the largest interaction of 
the DE mechanism, should induce a high spin ground state. However, the artificially large 
J makes that the hopping interaction is not large enough to restaure a ferromagnetic order 
in the usual DE model. As a consequence, the lowest potential energy well is the doublet 
state one. The minimum of this potential energy curve lies at q min =±0.97, 0.06 eV below 
the minimum in the octet curve at ±0.32. Owing to the antiferromagnetic contribution of 
the non-Hund states in the refined DE model, the optimized exchange integral J is properly 
estimated and both the spectrum at the experimental geometry and the potential energy 
curves of the different states reproduce the correct physics of the system. 
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FIG. 3: Potential energy curves of the eight states obtained from the optimized parameters of the 
usual DE model. The red line indicates the optimized value of q m j n in the octet ground state i.e. 
the geometry for which the ab initio spectrum has been calculated. The blue line indicates the 
position of the doublet state well. 

V. CONCLUSIONS 

The extraction of local electronic structure parameters in the CE-phase of the half-doped 
manganite Lao.sCao.sMnOs reveals that the local physics can accurately be described with 
a model Hamiltonian that includes the usual double exchange physics of mixed valence 
compounds plus an additional term due to the non-Hund states. The large t- value (1.3 
eV) indicates a rather mobile electron and hence small charge disproportionation. In fact, 
Mulliken population analysis associates a slightly higher charge to the Mn-ion in the Jahn- 
Teller distorted site. The difference between both sites is 0.07 electron. Although absolute 
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FIG. 4: Potential energy curves of the eight states obtained from the optimized parameters of the 
refined DE model. The red line indicates the optimized value of q m j n in the octet ground state 

Mulliken charges are not reliable, the relative values can be considered as indicative for 
tendencies. Moreover, the observed charge disproportionation of 0.07 is consistent with the 
maximum of 0.18 predicted in model Hamiltonian studies.— 

The PKS model provides a simple recipe to obtain the potential energy curves of the 
low-lying states as function of the distortion parameter q and hence obtain more detailed 
information about the local electronic physics. The curves for the low-lying states show a 
double well, where the minimum of the high spin ground state curve lies at q ± 0.38. This 
corresponds to a situation in between the completely delocalized solution with q min = 
and the completely localized solution for g m j n =±l. The latter solution has a dispropor- 
tionation of 1 electron while for the QW n =0 situation there is no barrier and the charge 
disproportionation has completely disappeared.— 
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The height of the barrier between left and right localized solutions is 0.04 eV for the 
octet ground state of the embedded cluster. Our model study only considers local distortions 
around the two Mn-sites considered in the cluster and the collective nature of the distortions 
in the CE-phase could lead to higher barriers. Obviously these effects should be included 



in order to be conclusive about the nature of the local electronic structure of the ha 



f- 



doped manganite; charge ordered or delocalized (Zener polaron) as proposed in Ref. |J 
Nevertheless, the smallness of the barrier is in line with the two different interpretations 
obtained by experimentalists and the controversy in the literature. 
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